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PROCEDE DE CALIBRAGE D^UN APPAREIL DMMAGERIE 
RADIOLOGIQUE NECESSITANT UN NOMBRE LIMITE D'ACQUISITIONS 

Linvention concerne le domaine medical et plus partlculierement le 
5 domaine des appareils d'imagerie radiologique. 

Ces appareils comprennent classiquement une source de rayons X 
et un detecteur d'images montes aux deux extremites d'un bras en C, entre 
lesquels est positionne Tobjet dont on souhaite realiser une image. Les 
images acquises par le detecteur sont le resultat d'une projection de Tobjet 
10 dans le plan du detecteur. 

Le bras en C est commande pour se deplacer autour de Tobjet et 
realiser ainsi plusieurs acquisitions correspondant a differents points de vue 
de I'objet. Des moyens de traitement permettent de reconstruire a partir des 
images bidimensionnelles acquises un modele tridimensionnelle de I'objet. 
15 Cette reconstruction suppose bien entendu de connaitre de maniere precise 
les differentes positions de Tappareil ainsi que ses caracteristiques 
geometriques. 

Le modele tridimensionnel obtenu peut etre utilise par le chirurgien 
avant une intervention pour se familiariser avec la partie de I'anatornie sur 

20 laquelle 11 doit pratiquer une intervention. II peut egalement etre utilise 
pendant Tintervention. A cet effet, le chirurgien dispose d'un equipement lui 
permettant d'afficher en temps reel des vues bidimensionnelles de la partie 
de Tanatomie, ces vues etant calculees a partir du modele tridimensionnel. 

La reconstruction du modele tridimensionnelle requiere un 

25 "calibrage geometrique" prealable de Tappareil d'imagerie. Ce calibrage 
permet de lier Tespace tridimensionnel aux informations bidimensionnelles 
fournies par les differentes projections bidimensionnelles. Si ce calibrage 
est grossier, la qualite du modele tridimensionnel reconstruit presentera des 
degradations. 

30 Une technique connue de calibrage consiste a disposer dans le 

champ des rayons X, des marqueurs positionnes sur un fantome qui. 
servent de reperes dans Tespace, et a realiser une serie d'acquisitions. Les 
positions de ces marqueurs dans Tespace tridimensionnel etant connues, 
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on peut en deduire la geometrie de Tacquisition pour chaque projection par 
inversion d'un systeme d'equations derive de la position des marqueurs sur 
les images projetees. Une telle technique est par exennple decrite dans le 
document US 5 442 674 (publie le 15 aout 1995). 
5 En general, lors de la realisation de Timage de la partie anatomique, 

Tappareil est commande pour realiser une s^rie d'acquisitions dans les 
memes conditions geometriques que la serie d'acquisitions de calibrage, de 
sorte la geometrie de chaque acquisition soit connue avec precision. Ceci 
est possible du fait que Ton consldere que les mouvements du bras en C 

1 0 sont repetables. 

Si toutefois les acquisitions de la partie anatomique a imager ne 
sont pas effectuees selon les memes conditions geometriques que les 
acquisitions de calibrage (nombre de vues, positions angulaires, Vitesse de 
rotation, positions initiales et finales), alors il est necessaire de realiser un 

15 nouveau calibrage, ce qui augmente le nombre de prises de vues realisees 
pour le calibrage. 

Par ailleurs, le nombre d'acquisitions necessaires pour obtenir un 
modele tridimensionnel de bonne quaiite depend du type de la structure 
anatomique a imager. Typiquement, la reconstruction tridimensionnelle de 

20 structures osseuses necessite environ 120 prises de vues tandis que la 
reconstruction tridimensionnelle de vaisseaux sanguins necessite environ 
40 prises de vue. 

II en resulte que le nombre d'acquisitions de calibrage doit toujours 
etre maximal de sorte que ce calibrage soit valable pour imager tout type de 

25 structures. 

Un but de invention est de minimiser le nombre d'acquisitions 
necessaires au calibrage de Tappareil. 

A cet effet, I'invention propose un procede de calibrage d'un 
appareil d'imagerie radiologique comprenant une source et un detecteur, 
30 comprenant les etapes consistant a : 

- deplacer I'appareil par rapport a un objet de calibrage et realiser 
une serie d'acquisitions, chaque acquisition etant associee a une position 
de calibrage de I'appareil, 
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- a partir des acquisitions realisees, determiner des parametres de 
projection associes a chaque position de calibrage de Tappareil, 

le precede etant caracterise en ce qu'il comprend en outre Tetape 
consistant a : 

5 . - pour una position supplementaire qui n*a pas ete prise par 

Tappareil lors de Tetape d'acquisition, determiner des valeurs de parametres 
de projection associes a cette position supplementaire en fonction de 
parametres associes aux positions de calibrage. 

Ce precede peut etre applique a un nombre quelconque de 
10 positions supplementaires, differentes des positions de calibrage. II a pour 
effet de multiplier par calcul le nombre de positions de calibrage obtenues. 

Ce procede permet par consequent de ne realiser qu'une quantite 
limitee d'acquisitions. 

En outre, ce procede evite de realiser des nouvelles acquisitions de 
15. calibrage dans le cas oCi les acquisitions realisees seraient en nombre 
insuffisant. 

Les parametres de projection peuvent etre des parametres 
geometriques caracteristiques du positionnement de rappareil ^dans 
Tespace (parametres extrinseques) ou des parametres geometriques 
20 caracteristiques de la source et du detecteur (parametres intrinseques). 

Ces parametres peuvent egalement etre des coefficients d'une 
matrice globale definissant la projection d'un objet de Tespace 
tridimensionnel dans le plan du detecteur. 

D'autres caracteristiques et avantages ressortiront encore de la 
25 description qui suit, laquelle est purement illustrative et non limitative et doit 
etre lue en regard des figures annexees parmi lesquelles : 

- la figure 1 represente schematiquement un appareil d'imagerie du 
type a bras en C, 

- la figure 2 represente schematiquement un bras en C supportant 
30 une source et un detecteur, 

- la figure 3 represente les differentes positions prises par la source 
lors des acquisitions de calibrage, 
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4 

- la figure 4 represente les positions supplementaires pouvant etre 
prises en compte pour I'application du procede de Tinvention, 

- la figure 5 est un bloc diagramme representant schematiquement 
les differentes etapes d'un mode de mise en ceuvre du procede de 

5 rinvention, 

- la figure 6 est un bloc diagramme representant schematiquement 
les differentes etapes d'une variante du mode de mise en ceuvre de la 
figure 5, 

- la figure 7 est un graphique illustrant le resultat d'un lissage des 
10 parametres de projection calcules par le procede de Tinvention. 

Sur la figure 1, Tappareil d'imagerie represente comprend de 
maniere classique une base 1 fixe par rapport au sol et sur laquelle est 
monte un bras en C 2. Le bras en C supporte a chacune de ses extremites 
une source 3 constituee par un tube a rayons X et un detecteur d'image 4. 
1 5 Ces elements sont disposes de sorte que le foyer de la source 3 et le plan 6 
du detecteur 4 se trouvent dans des positions diametralement opposees sur 
le bras en C, L'apparei! est positionne a proximite d'une table 5 sur laquelle 
sera installe un patient. Le plan de la table 5 s'etend entre les deux 
branches du bras en C, c'est a dire entre la source 3 et le detecteur 4. Le 
. 2.0 . bras en C 2 peut etre commande en rotation autour d'axes X ou Y, ces deux 
axes se coupant en un point O, centre de rotation du bras en C. 

Durant la phase de calibrage de Tappareil, un fantome comprenant 
des marqueurs est positionne sur la table 5 entre la source 3 et le detecteur 
4, sensiblement au niveau du point O. Durant cette phase, le bras en C est 
25 commande en rotation et Tappareil realise une serie de n acquisitions. 
Chaque image acquise est une projection des marqueurs, de I'espace 
tridimensionnel vers le plan 6 du detecteur 4. 

Sur la figure 2, on a represente schematiquement le bras en C 2 , le 
point focal S de la source 3 et le plan 6 du detecteur 4 ainsi qu'un objet 10 
30 dont on souhaite obtenir un modele tridimensionnel. 

Les parametres geometriques d'acquisition de Tappareil sont par 
exemple les suivants : 
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- six parametres extrinseques (trois translations, trois rotations) qui 
definissent la position dans l^espace du bras en C dans un repere (X, Y, Z) 
tridimensionnel, 

- cinq parametres intrinseques definissant ia projection dans le plan 
du detecteur : deux facteurs d'agrandissement (au. av) suivant deux axes 
d'un repere bidimensionnel (u, v) lie au plan, les coordonnee (uo, vo) de la 
projection s du point focal S de la source 3 ainsi qu'un angle de deviation 
que Ton considere nu! de manlere generate. 

Soil un point M de Tobjet a imager de coordonnee (Xm. Ym, Zm) dans 
Tespace et m son image projetee de coordonnees (Um, v^) dans le plan du 
detecteur. On a : 
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oCi R et T sont les matrices de rotation et de translation definissant 
la position du bras en C dans le repere tridimensionnel (X, Y. Z) et s est un 
facteur d'echelle. 

La matrice de projection P est done une matrice 3x4 dont les 12 
coefficients sont les resultats de combinaisons des 11 parametres (six 
parametres extrinseques et cinq parametres intrinseques) et d'un facteur 
d'echelle. 

Le diagramme de la figure 5 represente les differentes etapes d'un 
precede de reconstruction tridimensionnelle. Ce precede comprend une 
phase prealable de calibrage 100 et une phase d'acquisition 200 de la 
partie anatomique a imager. 

Selon une premiere etape 110 de la phase de calibrage, un fantome 
est positionne sur la table et le bras en C est entraTne en rotation. 

La figure 3 represente les differentes positions Xi, X2, X3, ...Xp prises 
par la source 3 autour du fantome 7 lors de la serie d'acquisitions de 
calibrage. 
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Selon une deuxieme etape 120 representee a la figure 5, des 
moyens de traitement calculent pour chaque position Xi, X2, X3, ...Xn de la 
source 3, une matrice de projection Pxi (definie par la relation [1] ci-dessus). 
La projection Pxj associee a une position X] de la source est determin6e a 
5 partir de la configuration des marqueurs sur Timage acquise. la position des 
nnarqueurs dans I'espace etant connue. 

Un algorithme de calibrage mis en oeuvre par les moyens de 
traitement permet d'associer a chaque position de calibrage Xi, X2, X3, ...Xn 
une projection Pxi, Px2, Pxa, ...Pxn. 
10 Selon une troisieme etape 130, les moyens de calcul determinent 

les parametres geom6triques axi, 3x2, . .axn d'acquisition de I'appareiL Pour 
chaque projection Pxi, ces parametres axi sont au nombre de 11 (six 
parametres extrinseques et cinq parametres intrinseques definis 
precedemment). 

15 Selon une quatrieme etape 140, les moyens de traitement estiment 

pour ces positions supplementaires yj des parametres geometriques 

d'acquisition ayi, ay2, ...ayn-i. 

Ainsi que represente sur la figure 4, on considere une serie de 

positions supplementaire yo, yi, y2, ...yn reparties sur la trajectoire de la 
20 source. La position yo est par exemple positionnee en amont de la position 

xi et la position yn est positionnee en aval de la position Xn. Chaque 

positions yj 0=1 . .n-1) est localisee sur la trajectoire de la source entre deux 

positions successives de calibrage x\ et X|+i (i<n). 

L'estimatiqn des parametres geometriques d'acquisition ayi, ay2, 
25 ...ayn-i est possible dans la mesure ou la trajectoire de la source est 

continue. II en resulte que la variation des parametres d'acquisition est 

egalement continue. 

A cet effet, les moyens de traitement realisent une interpolation des 

parametres geometriques axi, ax2, . axn determines pour les positions Xi, X2. 
30 X3, ...Xn. Chaque parametre geometrique ayj est done une combinaison des 

parametres axi - 
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n 

De maniere preferentielle, cette interpolation ne prend en compte 
qu'un nombre limite de positions de calibrage, inferieur a cinq. Par exemple, 
les moyens de traitement peuvent ne prendre en compte que les deux 
5 positions de la source Xj et Xj+i qui encadrent directement la position yj. 

Pour les positions extremes yo et yn, les moyens de traitement 
realisent une extrapolation des parametres geometriques axi, ax2, ...axn 
determines pour les positions xi, X2, xa, ...Xn. De la meme maniere, cette 
extrapolation ne prendra en compte qu'un nombre limite de positions de 
10 calibrage, inferieur a cinq. Par exemple, les moyens de traitement peuvent 
ne prendre en compte que les deux positions de la source Xj et Xi+i les plus 
proches de la position yj. 

Les interpolations et extrapolations realisees a partir de deux 
. positions peuvent avantageusement etre lineaires. '-"'^ 
15 Pour un nombre superieur de positions prises en compte, les 

fonctions d'interpolation peuvent etre des fonctions polynomiales, 
rationnelles ou d'un autre type. r 
Selon une cinquieme etape 150, les moyens de traitement 
deduisent de ces parametres les matrices Pyi, Pya, ...Pyn de projection 
20 associees aux positions supplementaire yo, yi, ya, -- .yn- 

La phase 100 de calibrage a pour effet d'augmenter le nombre de 
positions de calibrage obtenues par rapport au nombre n initiale de 
positions d'acquisitions. L'appareil d'imagerie est ensulte commande lors de 
la phase d'acquisition 200 pour realiser une image d'une structure 
25 anatomique d'un patient. 

Selon une premiere etape 210 de la phase d'acquisition, un patient 
est installe sur la table et le bras en C est entrame en rotation. Les 
acquisitions d'image peuvent etre realisees pour une partie ou pour 
Tensemble des positions de calibrage yo, Xi, yi, X2, yz, xa, ...Xn, yn de la 
30 source. 
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Selon una deuxieme etape 220, les moyens de traitement 
reconstruisent un modele tridimensionnel de structures anatomiques du 
patient a partir des images acquises et des projections Pxi, Pyi, Px2, Py2, 
Px3,.,.Pyn, Pxn associees a ces images. 
5 Le diagramme de la figure 6 represente une variante du procede de 

reconstruction tridimensionnelle de la figure 5. 

Le procede est similaire a celui de la figure 6 excepte que les 
coefficients Pyo, Pyi, Py2 ...Pyn des matrices Pyi, Pya, ...Pyn projection sont 
determinee directement par interpolation ou extrapolation des coefficients 
10 pxi, Px2, Px3 . .Pxn des matrices Pxi, Px2, Pxs, ...Pxn. 

Ainsi, a I'etape 150, les moyens de traitement realisent une 
interpolation des matrices Pxi, Px2, Px3, ...PXn determinees pour les 
positions Xi, X2, X3, ...Xn a Tetape 120. Chaque matrice Pyj est done une 
combinaison des matrices Pxj : 

15 Pyj=J]QjxPx, 

De la meme maniere que precedemment, cette interpolation peut 
ne prendre en compte qu'un nombre limite de positions de calibrage, 
inferieur a cinq. Par exemple, les moyens de traitement peuvent ne prendre 
en compte que les deux positions de la source Xj et Xi+i qui encadreht 

20 directement la position yj. 

Pour les positions extremes yo et yn, les moyens de traitement 
realisent une extrapolation des parametres geometriques axi. 3x2, . .axn 
determines pour les positions Xi, X2, X3, ...Xn. Cette extrapolation pourra ne 
prendre en compte qu'un nombre limite de positions de calibrage, inferieur a 

25 cinq. Par exemple. les moyens de traitement peuvent ne prendre en compte 
que les deux positions de la source x, et Xi+i les plus proches de la position 
Yi- 

Les interpolations et extrapolations realisees a partir de deux 
positions peuvent avantageusement etre lineaires. 
30 La variante de procede de reconstruction representee a la figure 6 

presente Tavantage de ne pas necessiter d'etape 130 de caicul des 
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parametres geometriques de projection, Tinterpolation etant realisee 
directement a partir des matrices de projection Pxi, Px2, PX3, ...Pxn 
determinees lors du calibrage. Par consequent, cette variante est simplifiee 
par rapport au procede de la figure 5 et conduit a des durees de traitement 
5 plus courtes. 

Toutefois, on notera que le procede represents a la figure 5 permet 
d'avoir acces aux parametres geometriques axi, 3x2, - axn. En particulier, ce 
procede autorise un lissage de ces parametres. 

Ainsi que represente en pointilles sur la figure 5, il est possible 
10 d'imposer des contraintes sur certains parametres prls en compte par 
Talgorithme de calibrage permettant de determiner les parametres 
geometriques des acquisitions de calibrage. 

On impose par exemple une loi definissant certains des parametres 
prls en compte par Talgorithme. Cette loi peut consister a fixer certains 
15 parametres ou a les definir par une fonction plus complexe dependant 
d'informations connues a priori et/ou d^autres parametres geometriques de 
projection. 

On impose par exemple que la position de la projection s de .la , 
source se trouve toujours au centre du detecteur, ce qui revient a fixer (uo, 

20 Vo). On peut egalement fixer la distance focale entre le point focal de la 
source, ce qui revient a fixer les facteurs d'agrandissement (au, av). 

On a pu constater que le fait de fixer Tun ou plusieurs de ces 
parametres permettait de regulariser les resultats obtenus et par 
consequent de recjuire le bruit qu'ils contiennent. 

25 La figure 7 illustre la variation d'un coefficient ty de la matrice T 

definissant la translation du bras en C selon Taxe Y pour un ensemble de 
positions definies par Tangle de rotation du bras en C, La courbe A 
represente les resultats obtenus lorsque aucun parametre geometrique 
n'est fixe. La courbe B represente les resultats obtenus lorsque les 

30 parametres (uo, Vq) ont ete fixes. 

Les precedes de reconstruction decrits precedemment permettent 
de calculer des parametres de projection pour des positions 
supplementaires de la source pour lesquelles il n'y a pas eu d'acquisition de 



1 er depot 



10 



calibrage. Ces positions supplementaires peuvent etre en nombre 
quelconque ou reparties de maniere quelconque sur la trajectoire de la 
source en fonction des besoins de la reconstruction tridimensionnelle finale. 
En particulier, il est possible d'intercaler autant de positions yj 
supplementaires entre les positions de calibrage que necessaire. 
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REVENDICATIONS 

1. Procede de calibrage d'un appareil d'imagerie radiologique 
comprenant une source (3) et un detecteur (4), comprenant les etapes 

5 consistant a : 

- deplacer I'appareil (2, 3, 4) par rapport a un objet de calibrage (7) 
et realiser une serie d'acquisitions, chaque acquisition etant associee a une 
position de calibrage (Xj) de Tappareil, 

- a partir des acquisitions realisees, determiner des parametres de 
10 projection (axi, Pxi) associes a chaque position de calibrage (Xj) de Tappareil, 

le procede etant caracterise en ce qu'il conriprend en outre I'etape 
consistant a : 

- pour une position supplementaire (yj) qui n'a pas ete prise par 
Tappareil lors de I'etape d'acquisition, determiner des valeurs de parametres 

15 de projection (ayj, pyj) associes a cette position supplementaire en fonction 
de parametres (axi, Pxi) associes aux positions de calibrage (Xi). 

2. Procede selon la revendication 1, caracterise en ce que les 
parametres de projection (axi, Pxi) comprennent des parametres 
geometriques caracteristiques du positionnement de rappareil dans 

20 Tespace. 

3. Procede selon Tune des revendications qui precedent, 
caracterise en ce que les parametres de projection (axi, Pxi) comprennent 
des parametres geometriques caracteristiques de la source (3) et du 
detecteur (4). 

25 4. Procede selon Tune des revendications qui precedent, 

caracterise en ce que les parametres de projection comprennent des 
coefficients (pxi) d'une matrice globale (PXj) definissant la projection d'un 
objet de Tespace tridimensionnel dans le plan (6) du detecteur (4). 

5. Procede selon Tune des revendications qui precedent, 

30 caracterise en ce que les parametres de projection (ayj, pyj) associes a la 
position supplementaire (yj) sont determines par une loi d'interpolation ou 
d'extrapolation des parametres de projection (axi, Pxi) associes a des 
positions de calibrage. 
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6. Precede selon la revendication 5, caracterise en ce que la loi 
d'interpolation est lineaire, polynomiale ou rationnelle. 

7. Precede selon Tune des revendications qui precedent, 
caracterise en ce que les parametres de projection (ayj, Pyj) associes a la 

5 position supplementaire (yj) sont determines par combinaison de 
parametres de projection (axi, Pxi) associes aux positions de calibrage (xi) de 
Tappareil les plus proches de la position supplementaire (Vj), le nombre de 
positions prises en compte etant inferieur a 5. 

8. Precede selon la revendication 7, caracterise en ce que le 
10 nombre de positions (Xj) prises en compte est egal a 2. 

9. Precede selon Tune des revendications qui precedent, 
caracterise en ce que lors de Tetape de determination des parametres de 
projection (axi, Pxi) associes a chaque position (Xi) de calibrage de Tappareil, 
certains parametres de projection sont definis par une loi dependant 

15 d'informations connues a priori et/ou d'autres parametres geometriques de 
projection. 

10. Precede de reconstruction d'un modele tridimensionnel a partir 
damages acquises par un appareil d'imagerie radiologique, caracterise en 
ce qu'il comprend une phase de calibrage (100) selon Tune des 

20 revendications 1 a 9 et une phase d'acquisition (200) comprenant les 
etapes consistant a : 

- deplacer I'appareil (2, 3, 4) par rapport a une structure (10) a 
imager et realiser une serie d'acquisitions, chaque acquisition etant 
associee a une position de calibrage ou a une position supplementaire (xj, 

25 Vj) de la source (3). 

- reconstruire un modele tridimensionnel de la structure (10) a partir 
des images acquises et des parametres de projection (axi, Pxi, ayj. Pyj) 
associes, determines lors de la phase de calibrage. 



30 
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